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For a system of spinless fermions in a disordered meso- 
scopic ring, interactions can give rise to an enhancement of 
the persistent current by orders of magnitude. The increase 
in the current is associated with a charge reorganization of 
the ground state. The interaction strength for which this 
reorganization takes place is sample-dependent and the log- 
averages over the ensemble are not representative. In this pa- 
per we demonstrate that the Hartree-Fock method closely re- 
produces results obtained by exact diagonalization. For spin- 
less fermions subject to a short-range Coulomb repulsion U 
we show that due to charge reorganization the derivative of 
the persistent current is a discontinuous function of U. Hav- 
ing established that the Hartree-Fock method works well in 
one dimension, we present corresponding results for persistent 
currents in two coupled chains. 
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In a normal-metal mesoscopic ring threaded by a mag- 
netic mix Q, measured values of the persistent current 
Q are two orders of magnitude larger than predicted. 
This result suggests that a quantitative theory must treat 
electron-electron interactions and disorder on an equal 
footing. Previous studies of spinless fermions concen- 
trated on the behaviour of ensemble averages of the per- 
sistent current and led to the conclusion that repulsive in- 
teractions cannot significally enhance the current. It was 
therefore argued that only systems which include spin 
could show such an increase of the current However, 
more recently it was shown Q] that for one-dimensional 
systems of spinless fermions interacting through a short- 
range Coulomb repulsion U, the persistent current does 
increase. This enhancement of the current is accompa- 
nied by a charge reorganization of the ground state which 
happens at different values of the interaction strength U, 
depending on the disorder realization, and therefore the 
ensemble averaged persistent current may not be rele- 
vant. 

The results of [Q] were obtained using the density ma- 
trix renormalization group (DMRG) technique Q which 
is essentially exact and contains correlation effects. How- 
ever, this technique is computationally demanding and 
not easily extended to higher dimensions. Therefore it is 



of interest to determine whether or not the results of Q] 
are contained in a mean-field description using a single 
Slater determinant. In this paper we address this ques- 
tion, by presenting results obtained using the Hartree- 
Fock method for spinless fermions interacting via a short- 
range Coulomb repulsion. In one dimension we show that 
the Hartree-Fock method agrees with the exact results of 
H , reproducing the predicted behaviour of the persistent 
current as well as the sample-dependent charge reorgani- 
zation. After establishing the validity of the method, we 
extend the calculation to a quasi-one dimensional system 
comprising two parallel chains. 

The total Hamiltonian for a system of N spinless 
fermions on a disordered chain (ID) of M sites is 
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The operators cj and c$ are creation and annihilation 
operators for an electron on site i, the on-site energies 
£j are random variables uniformly distributed over the 



interval 



•^f to and Uij is a nearest- neighbour in- 
teraction of the form 



Ui, 



U if j = % ± 1 
otherwise 



(2) 
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The hopping elements Uj are restricted to nearest 
neighbours with U^±i = t = 1 except at the ends of 
the chain, for which t±N = £/vi = 1 (periodic boundary 
conditions) or tiN = tjvi = —1 (anti-periodic boundary 
conditions). 

The single-particle Hartree-Fock equation correspond- 
ing to equation ([!]) is of the form 

- £ii_i* n (i - 1) - t ii+ i* n (i + 1) 

N M 

+ EEi* m wi 2 ^*"«- 

m— 1 j—1 

N M 

E E ** m 0')* m (0tf«* n cj) = £ n * n (o (3) 

m—1 j—1 

where \l/ n (i) is the amplitude of the nth single-particle 
wavefunction on site i. The third and fourth terms are 
the direct (Hartree) and exhange (Fock) potentials, re- 
spectively. 
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The parameters N, M, t, W, and U were chosen to be 
exactly the same as in Q so as to directly test the accu- 
racy of the method. Here we present results for M = 20 
sites and N = 10 particles (half-filling) and the strength 
of the disorder was taken to be large, W = 9. At zero U 
the charge density of a given sample is spatially inhomo- 
geneous due to the presence of such disorder. As U in- 
creases we observe a reorganization of the charge, and for 
large U, obtain a homogeneous configuration in which the 
particles are equally spaced. This reorganization happens 
at different values of the interaction strength for different 
samples. As an example, figure [l] shows the charge den- 
sity for a single sample in the free electron case (U = 0) 
and for U — 20, which is large enough to yield a periodic 
array of charges. 

In addition to the charge density we have studied the 
phase sensitivity D, which is a measure of the derealiza- 
tion effect mentioned above and is defined by 

D(U) = ^AE (4) 

Here AE = (-l) N (E g (0) - E g {n)) is the difference 
in the ground state energy between periodic and anti- 
periodic boundary conditions, where the Hartree-Fock 
ground state energy is 

N M 

E 9 = 2 + E ** n m^ n (j)} (5) 

n— 1 ij — 1 

with hij = EiSij + tij6ij±i. 

In agreement with we find peaks in log D at sample- 
dependent values of U, associated with reorganization of 
the ground-state charge density. For positive U, figure ^ 
shows the ensemble average of log D along with results for 
four individual samples. For negative U the mean field 
equations do not converge and no results were obtainable. 
Fig. [2] is in remarkable agreement with the exact results 
of Q. For example the average of log-D, exhibits a local 
maximum around U ~ t. Following Q we have also 
examined the relative increase of the phase sensitivity 
with respect to the free fermion case, r\ = log D(U) — 
logZ?(0), and in agreement with [Q, obtain a log-normal 
distribution for 77 (figure |3|) . 

Having established the validity of Hartree-Fock the- 
ory as a method for computing the charge density and 
phase sensitivity of one-dimensional rings, we now extend 
our analysis to two coupled one-dimensional chains. For 
a system consisting of two rings with 20 sites in each, 
and 20 spinless fermions in total, we again examine the 
case of strong disorder, W = 9. In figure ^ we show 
the charge density for the two rings. As we can see, the 
charge reorganization that was present in ID is also obvi- 
ous for two chains. For very strong U the particles local- 
ize in the odd and the even sites of the first and second 
chains respectively, as expected classically. Thus, one 



again obtains the derealization effect associated with the 
crossover from an Anderson to a Mott insulator. 

The phase sensitivity and the probability distribution 
of the relative increase r\ are shown in figures || and |(], 
which again demonstrate that interactions produce an in- 
crease in the fluctuations of the current. However, in con- 
trast with a single chain, the average of log D no longer 
possesses a local maximum and instead decreases mono- 
tonically with increasing U . 

In summary, the above results demonstrate that the 
effects discussed in [|j are contained in a single Slater- 
determinant ground state and are describable by mean 
field Hartree-Fock theory. By extending the analysis to 
two chains, we find that the maximum in the average 
of logD is no longer present, which suggests that this 
feature may be a peculiarity of strictly one-dimensional 
systems. The fact that Hartree-Fock theory is applicable 
in one dimension, where mean field theories are least ac- 
curate, indicates that in higher dimensions Hartree-Fock 
theory should be sufficient to describe the ground state of 
a system with electron-electron interactions and disorder, 
at least in the strong disorder limit. 
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FIG. 2. Phase sensitivity for half-filling (M = 20, N = 10) FIG. 4. Charge density for a sample with N = 20 particles 

in one dimension for four samples. The average has been in two rings ( 2 D) with M = 20 sites in each ring for free 
obtained for 500 disorder realizations fermions and for U = and U = 20 
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FIG. 5. Phase sensitivity for half-filling (M = 40, N = 20) 
in two dimensions for four samples. The average has been 
obtained from 400 disorder realizations 
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FIG. 6. Probability distribution of r\ = log D(2) - log D(0) 
in 2D calculated from 1000 samples. < r\ >= —1.59 and 
a 2 = 2.2 
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